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ABSTRACT 

We report on the weak lensing detection of a filament between two galaxy clusters at 
z = 0.55, CL0015.9+1609 and RX J0018.3+1618. We conduct weak lensing analysis of 
deep multi-band Subaru/Suprime-Cam images with Lensfit. The weak lensing sig¬ 
nals from the filament are contaminated by signals from the adjacent massive clus¬ 
ters and we statistically subtract the cluster component using two different methods. 
Both methods yield consistent shear profiles on the filament with > 2a significance 
and the average surface mass density of the filament is (E) = (3.20 ± 0.10) x 10 14 /t 
M.Mpc 2 , which is in broad agreement with previous studies. On-going surveys 
such as Hyper Suprime-Cam will identify more filaments, which will serve as a new 
probe of structure formation in the Universe. 

Key words: gravitational lensing: weak, large-scale structure of Universe 
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1 INTRODUCTION 


Recent observations have revealed that the Universe is 
dominated by unknown components c alled dark mat¬ 
ter a nd dark energy jHinsha w et all 120131: iRest et aD 
120141: IPlanck Collaboration et all 20151) . The theoretical 
model called A Cold Dark Matter (ACDM) model, which 
contains these components, has been confronted with 
severa l observational tests s uch as correla ti on fu n ctions 
(e.g., Blake & Bridlel 20051: iTegmark et all 120061; iHiitsil 
20061: iMarfn etal. 12013) and Minko wski functionals (e.g., 


Hikage et al 20081: Petri et ail 120131) . These statistical ob¬ 


servables are in excellent agreement with theoretical 
predictions from the ACDM model. However, these tests 
only use ensemble average properties of the matter dis¬ 
tribution. Complementary tests can be done with voids 
and filaments, which are important constitutes of the 
large-scale structure of the Universe. Voids are underdense 
regio ns whose properties hav e been studied theoretically 
(e.g.. iF urlanetto & Pirar J 120061: iKamionk owski et aD 120091; 

IHiguchi et all l20lT ~ Krause et alJ ~ 20131: Sutter et alT l2015h . 


* E-mail: yuichi.higuchi@nao.ac.jp 


On the other hand, filaments are slightly overdense regions 
connecting massive haloes. 

Large N-body simulations indicate that fil amen t s are 
form e d in the proces s of merging and acc retion dBond et all 
Il996t ISpringell l2005t IColberg et all l2005h . Therefore, fila¬ 
ments affect properties of nearly haloes. The properties of 
haloes residing with filaments are well studied through N- 
body simulations ' Co dis et all [20121; lAragon-Calvo & Yangl 
120141: IHiguchi et al 120141) . In order to find such filamen- 
tary structures, a variety of algorithms have been developed 

dMartinez & Saarll2002l;ISheth & Iainl2003l;ISousbie et alJ2006l: 

iMead et all2010h . 

Although we can see cosmic web in simulations, 
filamentary structures are hard to be detected ob serva- 
tionally due to their low overdensitv dSousbie et all I2OO8I; 
iHevmans et all 120081: IClampitt et aid l2014h . Some studies 


claimed the detection of filaments with X-rav tiBriel & Henrv 

1995 

;lTittlev & Henriksen 1200 HI Fane et alj 2002; Durret et al. 

2003 

; Werner 2008), although it is difficult to distinguish 


whether those signals come from filaments or from haloes 
at the edge of the filaments. 


While X-ray probes only hot intergalactic medium, weak 
gravitational lensing can probe not only baryon but also 
dark matter which occupies a dominant part of a filament 
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2 Y.Higuchi et al. 


mass. However, the detection of filaments with weak lens- 
ing have not been very successful. While weak lensing de- 
tectio ns of filaments b e tween clusters have been claimed 
(e.g., Clowe et al. 19981: iKaiser et~ai] 1 19981: iGrav et ahl l2002l: 
iDietrich et a 1.1 l2004h . it has also been suggested that these 
detectio n signals might have come from systema tics of the 
analysis dGavazzi et afll2004i:lHevmans et alJi2008 ~). Recently 
more secure detections of filaments with weak gravitational 
lensing have been r eported in the systems A bell 222/223 
dDietrich et alJ l2012h and MACSJ0717.5+3745 auzac et al] 
l2012h . However, a lager sample is needed in order to test 
the ACDM model with large-scale structures, and also to in¬ 
vestigate observational properties of filaments. In this paper, 
we report the new detection of a filament with weak gravita¬ 
tional lensing. 

This paper is organized as follows. In Section [2] we de¬ 
scribe the lensing analysis methods with a focus on the ba¬ 
sics of weak lensing and halo properties. In Section[3] we de¬ 
scribe our data analysis method, including the selection of 
background galaxies, the lensing analysis with Lens fit, the 
reconstruction of the filament mass distribution and the eval¬ 
uation of contaminations from clusters. In Section|4l we show 
the results of lensing analysis. We summarize our result in 
Section [5] 

In this paper, cosmological parameters are 

set accord in g to the_result of the Planck satellite 

dPlanck Collaboration et all l2013h : Hubble constant 
Hq = 67.3 km/s/Mpc, cosmological constant Q m = 0.315, 
dark energy Q A = 0.685. 


2 THEORIES FOR WEAK LENSING 

First, we summarize the theory of weak gravita¬ 
tional lensing. We follow the formalism outlined in 
iBartelmann & Schneider! d200ll) . 


2.1 Estimation of lensing signals from images 


Weak lensing signals are obtained from observed galaxy el- 
lipticities. However, convergence k( 0) and shear j( 9 ) are not 
directly related to galaxy ellipticities. The galaxy ellipticity of 
a lensed object depends on the reduced shear g defined as 


g(0) 


7(0) 

i - 40 )' 


(i) 


In the weak lensing limit, the reduced shear is related to the 
galaxy ellipticity £ as 

£ ~ £ s + g, (2) 


where e s is an intrinsic ellipticity of the galaxy. On the basis of 
the cosmological principle, we can assume that orientations 
of galaxies are random, i.e., the ensemble average of intrinsic 
ellipticities should be zero. Therefore, the average observed 
ellipticity is related to reduced shear as 

( £ ) = (g) ~ (7>/ (3) 

where (■ ■ ■} indicates ensemble average in a bin. 

We need to calculate the lensing depth of source galax¬ 
ies in order to estimate physical quantities such as cluster 
masses from weak lensing signals. The lensing amplitude is 


proportional to the distance ratio averaged over the popula¬ 
tion of source galaxies, 



dz 


dp Pis 

dz D s ' 


(4) 


where Di s and D s are the angular diameter distances from 
the lens to the source and from the observer to the source, re¬ 
spectively. dp/dz is the probability distribution function of 
redshifts of source galaxies. The critical projected mass den¬ 
sity is described with the lensing amplitude as 


E 


cr 


47tG 


,D, 



-l 


(5) 


where c is speed of light and G is gravitational constant. In 
this paper, we assume that all galaxies are located at this dis¬ 
tance. 


2.2 Model of filament profile 


The filament profile has been studied through both simu¬ 
lations and observation s dColberg et al.l 120051; Dietrich et alJ 
l2012l:lHlguchi et alj2014l) . These studies showed that the con¬ 
vergence profile of filaments is described as a function of dis¬ 
tance 9 from the line which connects between haloes on the 
sky plane (hereafter halo-halo axis) as dColberg et alj2005h 


k( 9 ) 


kq 

1 + (9/0 c ) 2 ' 


( 6 ) 


where k 0 describes the convergence value at 9 = 0 and 6 C 
is the scale length of a filament. For spherically symmetric 
objects, it is common to use reduced shear values defined 
relative to the centre of objects as 


( g+ \ = ( -cos2 cp —sin2 cp \ / gi \ (?) 

\ gx / v -sin2 <p coslip J \ g 2 ) ' 

where cp is the angle between axis-dj and 9. g + and g x are 
tangential and cross reduced shear components, respectively. 
However, filame nt s are m ore li ke axial symmetric. Therefore 
following iHiguchi et alJ j2014h . for filaments we define the 
"tangential shear" at each point as g + relative to the closest 
point on the halo-halo axis to the point. In this definition, 
the tangential shear profile of filaments takes negative values 
near the halo-halo axis and positive values at large radii (see 
Section [3.21 . 


2.3 NFW profile 

Weak lensing signals from filaments are contaminated by sig¬ 
nals from nearby clusters. We estimate the effect from clus¬ 
ters by using the cluster mass profile, referred as the Nav arro, 
Frenk and White (hereafter NFW) model l\’avarro et alJ 
Il997lk w hich has well been studied in simulations and obser - 
vations dGao et alj|2008l:l0kabe et ai1l2010l:l0guri et ai1l2012h . 
These studies indicate that the three-dimensional mass pro¬ 
file of clusters declines as r _1 inside the scale r s and de- dines 
as r~ 3 outside r s . The NFW profile is described as 


P(r) = 


r/r s (\+r/r s y ^ < rvir ) 

0 otherwise, 


( 8 ) 
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where r s and r v i r are the scale radius and virial radius, re¬ 
spectively. The characteristic density p s is defined as 


Ps = 


_Avir(z)p(z)4 ir _ 

3 [In (1 + c vir ) - c vir / (1 + c vir )[' 


(9) 


where A v ; r (z) is the non-linear overdensity predicted by the 
spherical collapse model. We parametrize the mass profile 
with two parameters. One is the concentration parameter: 


Cvir — 


and the other is the virial mass: 
M. 


4 71 

= -y A vir (z)io(z)r; 


,3 

vir* 


( 10 ) 


( 11 ) 


The ana lytical lensing profile of t he NFW profile is well 
studied (see [Wright & Brainerdl[2000il . In order to constrain 
the cluster density profiles, we fit the lensing signals with 
the two fitting parameters c v i r and M v ; r by minimizing a chi- 
square value defined as 


X 


2 


{g+,im(0;) g+,ana(6/)} 

l ^g+,im 


( 12 ) 


where and g+,ana (#;) are circular-averaged tangen¬ 

tial shear values at z'-th bin estimated from data and model, 
respectively, and c/g ( . is dispersion in z-th bin estimated 
from data. 



AR.A. [arcmin] 


3 DATA AND ANALYSIS METHOD 

In this section, we summarize the data, weak lensing analysis 
methods and the methods for subtracting the cluster compo¬ 
nents. 


Figure 1. Galaxy distributions at the cluster redshift. The filled points 
show galaxies in the photometric redshift range between 0.5 < 
Zphoto — & an d z photo + c < 0.6 where cr is the photo-z estimation 
error. The circles show the viral radii of two clusters estimated from 
profile fitting in Section[3] 


3.1 Data reduction 

We analyze the mass distribution around CL0015.9+1609 
(hereafter CL0016) and RX T0018.3+1618 (hereafter RX T0018) 
at redshift z = 0.55 ^Connolly et al.ll996l:lKodama et al.l2005l: 
iTanaka et al.ll2005l. l2009h . These previous studies show that 
the galaxy distribution between these two clusters exhibit fil¬ 
amentary structure (see Figure QJ. In the paper, we use the 
data takenjwith the Suba ru Prime Focus Camera (Suprime- 
Cam: lMivazaki et al]2002l) mounted on the Subaru Te lescope 
during the nights of 2 003 September 25-26 iKodama et alJ 
l2005l:lTfanaka et alJ2005h . We retrieve these data from the data 
archive system, SMOKA iBaba et alJl2002t) . The field of view 
of the Suprime-Cam is 34 x 27 arcmin 2 with the pixel scale of 
0.20 arcsec. It is observed through several broad-band filters 
( BVRi'z') down to limiting magnitude B = 26.9, V = 26.2, 
R = 26.0, i' = 25.9 and z' = 24.6, respectively. The seeing size 
is abo ut 0 ".65 (FWHM) . For m ore details, see lKodama et alJ 
d2005h and ITanaka et al J 120051) . 

The data are reduced with a modified version of the 
LSST stack, following the standard procedures of bias sub¬ 
traction, flat-fielding, distortion correction, sky subtraction 
and coadding to generate the stacked image from the mosaic 
images. The flat frames are constructed with the dithered sci¬ 
ence frames. The MAG-AU TO in SExtractor is used as a mea¬ 
sure of the total magnitude llBertin & Arnoutdl 19961) . 


For weak lensing analysis, a precise PSF model is cru¬ 
cial and we select PSF reference stars as follows. Objects 
are detected in the stacked z'-band image by using_SExtrac- 
tor. We use the same detection criteria used in iKodama et alJ 
( 20051 ). Stars are then selected in a standard way by identi¬ 
fying the stellar sequence in magnitude m,/ vs half light ra¬ 
dius r^ plane. To be specific, we apply 20.5 m[ < 23.5, 
r^ ^ 0.75 arcsec and more than 20c above the median sky. In 
addition, STAR-CLASS>0.95 and FLAG= 0 in SExtractor are 
also imposed. We set the fifth oder polynomial coefficients to 
fit the spatial variation of PSFs. For calculating the astromet¬ 
ric shifts between different shots, stars with 20c are used. 

Since contaminations by galaxies near the cluster red¬ 
shift dilute weak lensing signals, for weak lensing analysis 
we select only background galaxies using photometric red- 
shifts. We_use the_photometric redshift catalog constructed 
Tanaka et al. 1 2005 ) with the photometric redshift code 


Kodama et all Il999i) . Specifically, galaxies in photomet- 


m 
of 

ric redshift range Zph oto — c > 0.6 are selected as back¬ 
ground galaxies based on the photometric redshift catalogue, 
where c is an error for each photometric redshift estima¬ 
tion. The photometric red shifts of galaxies are compared 
with spectroscopical data i Hughes et all Il995l: iMunn et alJ 
Il997l: iHughes & Birkinshaw 19981) . The mean difference be¬ 
tween spectroscopic redshift z spec and photometric redshift 
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Figure 2. Photometric redshift distributions of galaxies. Solid and 
dotted lines show the distribution for all galaxies and background 
galaxies defined by z photo a > 0.6, respectively. 


is z photo ~ z S pec ~ 0.004. Our photometric redshifts tend 
to be under estimated for galaxi es with V — i' ~ 1.5 about 
Az ~ —0.1 (Kodama et al. fl999l l. In order to investigate the 
effect of this tendency, we estimate tensing signals with sev¬ 
eral different redshift selection criteria. However, we do not 

find any significant change of our result. _ 

We use the code Le nsfit ^Miller et al.ll2007l:lKitching et alJ 
120081: iMiller et al.l l2013h for estimating the lensing signals. 
Lensing signals are estimated with the i '-band mosaic data. 
We choose 40 x 40 pixel 2 as a postage stamp size in order to 
estimate galaxy ellipticities with Lensfit. To eliminate galax¬ 
ies that are not suitable for measuring lensing signals, e.g., 
blended galaxies, lensing signals are estimated with only 
galaxies with the flag fitclass = 0 in Lensfit. The total num¬ 
ber of galaxies used for weak lensing analysis is 28,_315. For 
each galaxy, we assign a weight defined by fe.e. jMiller et alJ 
l2013h 


£ 2 - la 2 

c max ^ u £ 


+ CTr 


pop 


(13) 


where cr t is the one-dimensional variance in ellipticity of the 
likelihood surface, Up 0 p is the one-dimensional variance of 
the distribution of ellipticity of the galaxy samples, and £ max 
is the maximum allowed ellipticity. The re duced shear in the 
n-th bin is calculated with the weig ht as JSeitz & Schneideil 
ll995l:lBartelmann & Schneiderll200lf) 


(§/') ( d n) 


-N, 


E=1 w g,i £ V,i 
E = 1 W gri 


(14) 


where pi indicates tangential or cross shear components. The 
convergence value at each point is obtained with galaxy el¬ 
lipticity as 

k( 0) = —E R e [#* (0 - 0 ( ) w g/i £i\ /Yu w g,i' ( 15 ) 

lip- 7T L J 

& l l 


where n g is a number density of galaxies, and & is the kernel 


defined as 



( 16 ) 


where 0 S is a smoothing scale. The statistical error in the 
shear measurement for each bin is computed from the 
weighted averag e of the variance of the shear a 2 (see 
lOkabe et alJ I201C : lOguri et~aD l2012h . For this purpose, we 
identify 20 neighboring galaxies in the magnitude-galaxy 
scale length plane. The variance, which is defined as cr 2 = 
Ugj + a 2 2 , is computed over the neighboring samples. Using 
the statistical error of shear for the i-th galaxy a g i, the uncer¬ 
tainty on the tangential shear in each bin is estimated with 
the weight as 


cr g+^") 


lX^I^ 

2 fawj ) 2 


(17) 


The factor 1/2 accounts for the fact that tr ■ is the rms 
for sum of two distortion components. We assume that the 
correlation between different bins is negligibly small, i.e., the 
main source of the statistical error comes from the intrinsic 
ellipticities. 

In order to quantify the significance of the filament de¬ 
tection, we estimate the difference of the profiles from a null 
profile, i.e., g+,ana (0;) = 0 in equation (T2). We estimate the 
significance by comparing the estimated chi-square values 
with the chi-square distribution. In this paper, we use fila¬ 
ment shear profiles up to 15 arcmin. 

In the shear measurement, we select galaxies in the pho¬ 
tometric redshift range z — a ^ 0.6 and estimate the lens¬ 
ing depth for the galaxies with equation 0. The value of the 
lensing depth (Di s /D s ) is 0.419. Figure |2] shows the redshift 
distribution of galaxies. 


3.2 Estimation of the effect of the clusters 

First, we select the Brightest Cluster Galaxy for each cluster 
(hereafter BCG) and define a line connecting both BCG po¬ 
sitions, which we call BCG-BCG axis. Then, we define a fil¬ 
ament region, which are the region between both BCGs and 
outside a radius R; from each BCG (shaded region in Fig- 
ure|3]i. The filament length l is defined on the sky plane as 

^2d = E r ; + z ' ( 18 ) 

i 

where L 2 d is the distance between both BCG positions and 
the index labels the clusters. As explained in Section 12.21 
shear profiles of the filament region is calculated as a func¬ 
tion of distance from the BCG-BCG axis. 

Since lensing signals from filaments are weak, clusters 
at both edges of the filament can affect lensing signals of the 
filament significantly. Therefore, it is important to subtract 
such cluster contributions. In order to subtract the effect, we 
estimate a background which comes from cluster weak lens¬ 
ing signals, and subtract the background from the lensing 
signals in the filament region. The background is estimated 
with two different methods to check robustness against back¬ 
ground subtraction methods. The first method uses analyti¬ 
cal profiles of the clusters. The second method estimates the 
background from the opposite sides of the filaments. Figure[3] 
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Figure 4. Left panels show the Subaru/Suprime-Cam image and contours of the convergence field for CL0016 (upper) and RX J0018 (lower). 
The contours (solid line) show the significance above the mean of the field edge spaced in a unit of 0.5(7 from 0.5(7 up to 3(7 for clarity. The dashed 
circles indicate the viral radii estimated from the shear profile fitting in Section[3] Right panels show the radial profiles of tangential and cross 
shear as a function of distance from the each BCG position. The error bars show the 1(7 measurement errors. The solid lines show the best-fit 
NFW profiles. In the tangential shear profile for RX J0018, the values of the fifth, seventh and eighth bins are not displayed due to low values. 


illustrate the procedure which we explained in this subsec¬ 
tion. 


3.2.1 Subtracting the analytical profiles 

The first background subtraction method is to subtract clus¬ 
ter lensing signals estimated from the best-fit analytical 
model. Assuming that the lensing signals of the clusters are 
well modeled with the spherically symmetric NFW profile, 
the effects from both clusters to the lensing signals of the fil¬ 
ament region can be removed by subtracting the lensing sig¬ 
nals constructed with the best-fit NFW models. 


We fit tangential shear profiles of both clusters with 
the NFW profile by minimizing chi-square values defined 
in equation i fhZt . Figure [4] shows the tangential shear pro¬ 
files and the best-fit NFW profiles. The best-fit parameters 
are M vir = (2.09+““) x lO^/i^M© and c vir = 5.56+““ 

for CL0016, and M vir = (4.51+““) x 10 14 /r “M® and 

c v ir = 4.61+j'jg for RX J0018. The best-fit masses are con¬ 
sistent wihi_die_jesultsobtained_from_hie2Gray observa¬ 
tions faughes et alJ Il995l; ISolovveva et alJ 120071) . The two- 
dimensional lensing shear maps of the clusters are con¬ 
structed with these parameters. Then, for each background 
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filament 


background region background 



Figure 3. Definition of the regions and explanation for the back¬ 
ground estimation. The black points show the locations of BCGs. A 
region between clusters and outside R from each BCG are defined as 
the filament region (shaded region). In the first method for a back¬ 
ground estimation, shear profiles inside 3R V ir from BCGs (dotted re¬ 
gions) are estimated with the analytical models and subtracted. In 
the second method, the regions on the side opposite to the filament, 
where regions inside R from the each BCG are excluded, are defined 
as background (striped regions). Ellipses in the filament region indi¬ 
cate the distorted directions of images of background galaxies by a 
mass distribution of the filament. 


galaxies inside 3R v i r from BCG positions shear values of 
the cluster components are subtracted to remove the cluster 
contribution. We also estimate the cluster masses only with 
galaxies which do not exist in the filament region. The esti¬ 
mated parameters are almost consistent within 1 a. 


3.2.2 Empirical background subtraction using opposite sides of 
the clusters 

The second background estimation method is to subtract the 
the shear profile on the opposite sides of the clusters. An ad¬ 
vantage of this second approach is that it takes account of 
elongations of clusters which are ignored in the first method. 
As shown in Figure [3] (stripe region), the background is set 
to the same length as that of the filament length. The regions 
inside R, from each BCG are excluded. Then, the tangential 
shear in the background is estimated as a function of distance 
from the line along the BCG-BCG axis, which is subtracted 
from the tangential shear profile in the filament region to ob¬ 
tain the final profile. 


4 RESULT 

Figure [5] shows the convergence map. We find the mass 
bridge between CL0016 and RX J0018, which is also seen in 
the galaxy distribution (FigureQ]. 

Figure [6] shows the tangential and cross shear profiles of 
the filament for both background subtraction methods. The 
shear profiles are shown for different definitions of the ex¬ 
cluded region (R = 1 /4r v ; r , l/2r vi r and 3/4r v i r ) to check the 
sensitivity of our results on the size of the excluded region. 



AR.A. [arcmin] 

Figure 5. The mass map reconstructed with weak lensing. The clus¬ 
ter near the centre is CL0016. RX J0018 is located on the lower-right 
corner. The contours show significance above the mean of the field 
edge spaced in a unit of 0.5(7 from 0.5(7. 


We find that the tangential shear profiles are similar to the 
shear profile exp ected from t he conv ergence profile in equa¬ 
tion © (see also iHiguchi et alJl201~3) . whereas B-mode pro¬ 
files are consistent with zero for both background subtrac¬ 
tion methods. Table Q] summarizes the values of the signifi¬ 
cances for the tangential shear profiles up to 15 arcmin for 
different definitions of the excluded region, estimated from 
X 2 for the null shear profile. We find that the filament pro¬ 
files are detected at the significance > 2a for the case of 
empirical background subtraction method. For the analyt¬ 
ical background subtraction method, the detection signifi¬ 
cances are much higher, but B-mode profiles in the inner re¬ 
gion slightly deviates from zero, which may indicate insuf¬ 
ficient background subtraction due to the aspherical profiles 
of the clusters. To illustrate this point, we compute the de¬ 
tection significances also for the B-mode (cross) shear pro¬ 
files, which are summarized in Table [2] While the x 1 values 
are consistent with those expected from the x 1 distribution 
for the empirical background subtraction method, x 2 values 
are significantly larger than expected for the analytical back¬ 
ground subtraction method, suggesting significant residuals 
for the analytical background subtraction. 

We also check the dependence of the profiles on the size 
of bins and the selection of galaxies. For instance, we change 
the size of bins and change the photometric cut to Zph oto — 
a > 0.7 or Zph oto — cr ^ 0.8. However, these changes do not 
strongly affect our result. 

We also estimate the average projected mass of the fila- 
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Figure 6. Tangential profiles of the filament for different excluded regions (upper left: l/4r V ir; upper right: l/2r v j r ; lower: 3/4r V ij). The crosses 
show the values estimated with the analytical background subtraction method. The squares show the values for the empirical background 
subtraction method. The error bars show lot The solid and dotted lines show the best-fit profiles with 0 C = 2 arcmin for the analytical and 
empirical background subtracted profiles. 


Table 1. Significance for the tangential shear profiles up to 15 arcmin for different definitions of the excluded region. Column (1): background 
estimation method; Column (2)-(4): value of significance. 


background subtraction method 

l/4r vir 

size of the excluded region 
l/2r vir 

3/4r vir 

analytical model 

13.5c- 

7.36 a 

4.910- 

empirical 

5.95 a 

2.05 a 

1.75o- 


Table 2. Similar to Table[l] But for the cross shear profiles. Column (1): background estimation method; Column (2)-(4): value of significance. 


background subtraction method 

l/4r vir 

size of the excluded region 
l/2r vir 

3/4r vir 

analytical model 

8.29o- 

8.10 a 

3.30o 

empirical 

l.llo- 

0.25 a 

0.65 a 


© 0000 RAS, MNRAS 000, 000-000 







































8 Y.Higuchi et al. 


Table 3. Average projected mass density of the filament. The mass density is estimated with convergence values within 3 arcmin from the 
BCG-BCG line for each definition of the excluded region. The errors indicate lcr error. 


Definition of the excluded region 

l/4r vil 

l/2r vir 

3/4r vir 

Projected mass density [!O 14 /!M 0 /Mpc 2 ] 

3.23 ± 0.10 

3.21 ±0.10 

3.20 ±0.10 


Table 4. Best-fit parameters with 9 C = 2 arcmin. Column (1): background estimation method; Column (2)-(4): value of kq. 


background subtraction method 

l/4r vir 

size of the excluded region 
l/2r vir 

3/4r vir 

analytical model 

0.33 

0.29 

0.33 

empirical 

0.28 

0.19 

0.27 


ment from the convergence map, using convergence values 
up to 0o = 3 arcmin from the BCG-BCG axis. Here we use the 
convergence map before the background subtraction. Table[3] 
shows the average and error of projected mass density for 
different excluded regions. We find that the value of the pro¬ 
jected mass density is almos t consistent with previous weak 
tensing studies of filaments dDietrich et alj|20l3:llauzac et alJ 
12012b . We fit the tangential shear profiles with the analytical 
profile defined with equation (|6). For simplicity, we fix 9 C = 2 
arcmin. Table[4]shows the best-fit values for kq. Table[5]shows 
the projected mass within 3 arcmin from the BCG-BCG axis, 
which are estimated by equation (|6) with the best-fit param¬ 
eters. The estimated mass densities are higher than those es¬ 
timated from the convergence values that are estimated from 
the smoothed mass map. The fitting with the fixed parameter 
9 C might also be a part of reasons for this discrepancy. 


5 CONCLUSION 

In this paper, we have investigated the weak tensing prop¬ 
erties of a filament between the cluster CL0015.9+1609 and 
RX J0018.3+1618 with the Subaru/Suprime-Cam data. We 
have selected background galaxies with photometric redshift 
and estimated the tensing signals with Lens fit. We have sub¬ 
tracted the effects of the clusters from the tensing profile of 
the filament with two different methods. One is to subtract 
the best-fit NFW profiles, and the other is to subtract the pro¬ 
file obtained on the opposite sides of the filament. 

We have fitted the tangential shear profiles of the 
clusters with the NFW model, and obtained M V u. = 
(2.09+°'“) x IO^/i^Mq and c vir = 5.56+°;“ for CL0016, 

and M vir = (4.51++“) x 10 14 /i 'Mq and c vir = 4.61±]J° 
for RX J0018. The estimated viral radii are used to define the 
size of the excluded regions. 

We have found the mass bridge between two clusters on 
the convergence map. We have estimated the average surface 
mass density for three different excluded region definitions. 
The surface mass density for the excluded region definition 
of R = 3/4r vir is found to (3.20 ± 0.10) x lO 14 //M 0 Mp C - 2 
from the convergence map. This is consistent with previous 
studies of filaments. 


We have derived the shear profiles of the filament using 
different excluded region definitions and background sub¬ 
traction methods. For all the definitions of the excluded re¬ 
gion, we have found that the cross shear values are almost 
consistent with zero and tangential shear profiles have same 
profiles as those expected from simulations, indicating that 
the tensing signals of the filament are detected. We have com¬ 
pared the shear profiles with a null profile up to 15 arcmin 
from the BCG-BCG axis to find that the tensing signals of the 
filament are detected at the significance > lcr in all cases. 

The detection of filaments with weak gravitational tens¬ 
ing is difficult due to small values of tensing signals. How¬ 
ever, large-scale surveys such as Hyper Suprime-Cam and 
Large Synoptic Survey Telescope may enable direct observa¬ 
tions of more filaments. Statistics of weak tensing properties 
of filaments therefor e serve as a ne w tes t of structure forma¬ 
tion in the Universe jHiguchi et al.l[2014h . 
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